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ABSTRACT 


Extensive  steady-state  and  dynamic  design  information  has  been 
generated  for  tilting  pad  journal  bearings  with  emphasis  on  miniature 
cryogenic  turbomachinery  applications.  Coverage  includes  bearing 
length  to  diameter  ratios  of  from  0.5  to  2.0  and  values  of  the  com¬ 
pressibility  parameter  from  1.5  to  20.  Information  was  put  into  a 
format  for  use  by  design  engineers  and  a  separate  section  discusses 
design  procedures  plus  practical  considerations.  An  experimental 
program  was  completed,  using  room  temperature  air  as  the  lubricant, 
which  substantiated  the  analytical  information. 

Variable  grid  computerized  analysis  was  required  to  avoid  numerical 
problems  at  high  values  of  the  compressibility  parameter  A.  These 
methods  proved  very  successful  and  represent  a  significant  improvement 
over  constant  grid  techniques. 
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1 .  INTRODUCTION 


1.1  Primary  Objective 

The  whirl-free  characteristics  of  tilting-pad  journal  bearings  es¬ 
tablished  them  as  the  favored  candidates  for  use  in  miniature  cryogenic 
turbomachinery.  (See  Reference  1.)  A  limited  amount  of  design  informa¬ 
tion  for  this  bearing  was  developed  in  a  previous  program  described  in 
Reference  1.  The  objective  of  the  present  effort  was  to  produce  extended 
verified  design  information  for  the  tilting-pad  journal  bearing  over  a 
sufficient  range  to  cover  operating  conditions  anticipated  for  future 
hardware . 

1.2  Basic  Configuration 

The  tilting-pad  journal  bearing  is  schematically  shown  on  Figure  1. 

It  is  a  three-pad  bearing,  and  the  pad  pivots  are  equally  spaced  around 
the  pivot  circle.  The  tri-pad  arrangement  ensures  that  pivot  positions 
remain  on  a  single  circle.  One  pad  is  generally  spring-loaded  to  main¬ 
tain  a  specified  pre-load  and  allow  for  centrifugal  and  thermal  expansions 
of  the  shaft.  The  information  developed  treats  both  the  fixed  pivoted 
bearing  and  the  spring-loaded  configuration.  The  position  of  the  pivot  was 
fixed  at  65%  from  the  leading  edge  of  the  pad,  since  this  position  is  optimum 
with  respect  to  load  capacity  over  a  wide  operating  range.  Also,  the  over¬ 
all  pad  angle  was  100  degrees  for  all  investigations. 

1.3  Range  of  Data  Generated 

The  original  study  (Reference  1)  covered  values  of  the  compressibility 
parameter  (A)  up  to  3.5.  Because  of  possible  variations  in  the  clearance/ 
radius  ratio,  speed,  ambient  pressure  and  lubricant  viscosity,  and  also  to 
design  away  from  unstable  regimes,  it  was  decided  to  extend  the  range  of 
A  to  20.  Performance  for  additional  length  to  diameter  ratios  from  that  of 


Figure  L  Schematic  Representation  of  Pivoted-Pad  Journal  Bearing 
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Reference  1  was  also  accomplished.  Figure  2  is  a  summary  chart  that  shows 
the  range  of  variable  covered  for  the  steady-state  analysis  and  Figure  3 
depicts  the  parameters  covered  by  the  dynamic  analyses.  It  was  not  orig¬ 
inally  anticipated  to  do  dynamic  analyses  for  A  =  5,  but  examination  of  the 
data  between  A  =  3.5  and  A  =  10  indicated  that  an  interim  value  would  be 
beneficial . 

1.4  txperimental  Program 

The  primary  purpose  of  the  experimental  program  was  to  verify  analyt¬ 
ical  predictions,  both  steady-state  and  dynamic.  This  was  accomplished 
on  a  comparatively  large  size  rig  to  facilitate  installation  of  instru¬ 
mentation  and  measurement,  using  air  as  the  gaseous  medium. 

1.5  Design  Methods 

A  secondary  objective  of  the  program  was  to  present  the  data  in  a 
format  that  could  readily  be  used  by  design  engineers.  A  separate  section 
(Section  4)  is  devoted  entirely  to  design  procedures  and  practical  con¬ 
siderations,  for  the  benefit  of  design  engineers  charged  with  the  respon¬ 
sibility  of  establishing  gas-lubricated  journal  bearing  configurations  for 
cryogenic  applications. 
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KEY :  R  =  Shaft  radius 

L  =  Bearing  length 

D  =  Shaft  diameter 

SP  =  Single  pad 

FB  =  Full  bearing  (3  pads) 

/  =  Accomplished  under  original  contract 
X  =  Accomplished  under  present  contract 

Full  bearing  data  generated  for  displacements  directly  between  pads  and 
directly  at  pads  for  both  fixed  pivot  configuration  and  for  single  spring 
loaded  shoe. 


Figure  2.  Steady-State  Summary  Chart 
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KEY :  R  *  Radius  of  shaft. 

L  «  Bearing  length. 

D  *  Shaft  diameter. 

H  *  Dimensionless  pivot  film  thickness. 
p  *  h  /C  vher*-.  h  »  pivot  film  thickness  and 
Cp*  radial  clearance. 

/  ■  Accomplished  under  initial  contract. 

X  »  Accomplished  under  present,  contract. 


Figure  3.  Dynamic  Summary  Chart 


(a  =  100°;  <p/a  =  0.65) 
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2.  ANALYSIS 


2.1  Steady-State  Analysis 

The  analytical  procedure  for  obtaining  steady-state  performance 
is  comprehensively  described  in  Reference  1.  For  completeness)  portions 
are  repeated  and  elaborated  upon  in  Appendix  I.  Since  a  major  objective 
of  this  program  was  the  generation  of  design  data,  it  was  decided  to  auto¬ 
mate  the  analytical  procedures  so  that  a  large  amount  of  performance 
information  could  be  economically  produced.  The  sequential  steps  were 
as  follows: 

2.1.1  Production  of  Single-Pad  Steady-State  Data 

The  foundation  of  the  steady-state  analysis  is  determining  the  re¬ 
lationship  between  the  pivot  film  thickness  and  performance  parameters 
such  as  load  coefficient,  viscous  friction  coefficient  and  minimum  film 
thickness  for  a  single  pad.  This  data  is  then  transposed  into  full  three 
pad  bearing  performance.  Appendix  I  describes  the  computerized  analyses 
used  to  generate  the  single-pad  data  for  a  known  bearing  size  and  pivot 
position. 

2.1.2  Omp1£fion  of  Curve  Fitting  Routines  for  Single-Pad  Steady-State 
Data 

After  single  pad  data  is  generated  as  a  function  of  the  pivot  film 
thickness,  it  is  necessary  to  put  results  in  a  convenient  form  for  further 
use.  By  utilizing  curve-fitting  routines,  it  was  possible  to  establish 
polynomial  expressions  over  varying  ranges  of  the  pivot  film  thickness 
that  accurately  described  the  data  from  the  computer  program  used  to 
establish  single-pad  performance. 
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2.1.3  Generation  of  Computer  Program  to  Determine  Full-Bearing  Performance 

After  single-pad  data  was  described  by  analytical  expressions  as  a 
function  of  pivot-film  thickness,  as  described  in  paragraph  2.1.2  above, 
it  was  not  difficult  to  assemble  another  program  whose  function  was  to 
produce  full  bearing  data.  A  total  of  four  conditions  were  considered, 
as  follows: 

1.  Three-pad  bearing,  fixed  pivots,  displacement 
between  shoes. 

2.  Three-pad  bearing,  fixed  pivots,  displacement 
directly  at  pivot  of  one  pad. 

3.  Three-pad  bearing,  one  pad  spring  loaded,  dis¬ 
placement  directly  between  two  fixed  pads. 

4.  Three-pad  bearing,  one  pad  spring  loaded,  dis¬ 
placement  directly  at  pivot  of  spring  loaded 
pad. 

Performance  was  established  as  a  function  of  eccentricity  ratio 
and  pivot  circle  pre-load.  Appendix  I  contains  the  analytical  procedure 
for  converting  from  single  pad  to  full  bearing  performance  and  also  contains 
the  complete  set  of  steady-state  performance  curves.  The  use  of  these 
curves  is  described  in  Section  4,  Design  Procedures. 

2.1.4  Automatic  Plotting  of  Full  Bearing  Performance 

An  automatic  plotting  routine,  which  generated  performance  curves 
for  the  full  bearing  proved  successful  and  was  extensively  applied.  The 
computer  plots,  accomplished  by  a  Calcomp  X-Y  plotter  allowed  great 
savings  in  manpower  and  time  in  production  of  the  design  curves. 

2.2  Samples  of  Steady-State  Results 

Figure  4  shows  single  pad  load  coefficient  Vs.  pivot  film  thickness 
for  values  of  the  compressibility  ratio  A,  of  1.5  and  3.5  It  is  essen¬ 
tially  this  data  that  is  converted  to  polynomial  format  and  coefficients 
determined  by  insuring  a  very  accurate  curve  fit.  Then,  by  utilizing 
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Figure  4.  Individual  Pad  Load  vs.  Pivot  Film  Thickness 
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procedures  described  in  Appendix  I  full  bearing  data  is  determined  as 
a  function  of  eccentricity  ratio  e'  and  the  pad  pre-load  parameter  C'/C. 
From  this  data  curves  are  automatically  generated  (from  additional  curve 
fitting  routines)  to  produce  performance  plots  as  shown  on  Figure  5.  Here 
the  full  bearing  load  coefficient  is  plotted  as  a  function  of  bearing 
eccentricity  ratio  for  various  pre-loads.  Thus,  for  a  given  design  load, 
proper  pre-load  to  satisfy  a  specific  eccentricity  (which  can  be  related 
to  pivot  and  minimum  film  thickness)  can  be  established.  The  curves  also 
give  pivdt  film  thickness  on  all  three  shoes.  The  series  of  curves  shown 
on  Figure  5  also  permit  determining  off-design  performance.  Figures  6  and 
7  show  the  ether  two  parameters  of  primary  interest  namely,  pad  friction 
coefficient  and  minimum  film  thickness.  This  information  was  developed 
for  individual  pads  and  not  the  whole  bearing.  For  friction,  it  is  neces¬ 
sary  to  accumulate  the  effect  of  each  separate  pad. 

There  are  quite  a  voluminous  number  of  curves  similar  to  Figures  4 
through  7  covering  the  range  and  conditions  described  in  Section  1. 

Their  use  for  determining  steady  state  performance  will  be  shown  by  example 
in  Section  4,  Design  Procedures. 

2.3  High  Lambda  (a)  Analysis 

Numerical  difficulties  were  encountered  because  of  the  rapid  pres¬ 
sure  changes  in  the  bearing  trailing  edge  film  at  high  A's.  These  changes 
were  so  severe  that  dividing  the  bearing  circumference  with  a  uniform  grid 
spacing  resulted  in  the  computer  solution  becoming  unstable  or  converging 
to  the  wrong  solution.  By  designing  a  system  with  variable  grid  spacing, 
this  convergence  problem  was  eliminated.  The  technique  merely  increases 
grid  spacing  where  the  pressure  changes  slowly  (i.e.,  leading  edge)  and 
decreases  grid  spacing  where  the  pressure  changes  rapidly  (i.e.,  trailing 
edge) . 

The  steady-state  portion  of  the  computer  program  was  modified  to 
incorporate  this  variable  grid  technique.  Figure  8  is  a  plot  of  the  non- 
dimensional  pressure  distribution  in  the  circumferential  direction, 
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PIVOT  FILM  THICKNESS 


PAD  FRICTION  COEFFICIENT,  F-  Ff/PaCL 


VARIABLE  GRID  (0.18  TO  0.01) 


Figure  8.  Circumferential  Pressure  Distribution,  Variable  vs.  Fixed  G 


(i.e.,  in  the  direction  of  rotation)  taken  at  the  center  of  the  pad.  A 
small  lambda  value  (A  =  1.5)  was  used  to  check  the  validity  of  the  variable 
grid  approach.  The  solid  line  is  the  pressure  distribution  using  a  fixed 
grid  (circumferential  increment  *  0.0727  radians).  The  values  of  pressure 
signified  by  X  and  0  are  values  using  two  different  variable  grid  spacings. 

0  range  from  0.18  radians  at  leading  edge  to 
0.01  radians  at  trailing  edge. 

X  range  from  0.10  radians  at  leading  edge  to 
0.02  radians  at  trailing  edge. 

A  high  A  test  case  was  run  to  check  out  the  modified  variable  grid 
technique.  Figure  9  is  the  pressure  distribution  at  A  *=  20  using  fixed 
and  variable  grid  spacings.  The  dotted  line  uses  a  fixed  grid  and  the 
solid  line  uses  a  variable  grid.  The  variable  grid  generates  the  correct 
pressure  distribution  while  the  fixed  grid  generates  a  pressure  profile 
obviously  suffering  from  numerical  instability.  Figure  10  presents  the 
corresponding  film  thickness  distribution.  Appendix  II  describes  the 
theory  for  implementation  of  the  variable  grid  concept. 

2.4  Dynamic  Analysis 

The  big  advantage  of  tilting-pad  gas  bearing  as  compared  to  other 
types  is  their  relative  freedom  from  self-excited  fractional  frequency 
whirl.  However,  for  the  unusually  severe  conditions  encountered  with 
miniature  cryogenic  turbomachinery,  shoe  flutter  could  cause  problems. 
Reference  2  describes  an  analytical  technique  that  determines  whether 
flutter  will  be  excited  and  cause  bearing  failure.  The  technique  is 
called  the  "Step-Jump"  method.  The  bearing  is  assumed  to  be  in  its 
normal  equilibrium  position  as  if  it  were  to  be  stable.  The  journal  and 
pads  are  then  given  sequential  slight  disturbances  in  their  various  de¬ 
grees  of  freedom.  Since  the  fluid  film  is  compressible  the  response  to 
the  disturbances  is  not  immediate,  but  occurs  over  some  discrete  time 
interval.  The  responses  are  obtained  by  numerical  solution  of  the  time- 
dependent  compressible  Reynolds'  equation;  thus,  they  are  obtained  very 
accurately.  They  are  then  transposed  to  a  polynomial  form  so  they  can 
be  continuously  re-used  in  the  equations  governing  the  dynamics  of  the 
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system.  Then,  inertias,  masses  and  spring  properties  of  the  dynamic 
system  can  be  changed  ar.d  the  equations  solved  without  resolving  Reynolds' 
equation  for  the  fluid  film  iorces.  It  is  the  solution  of  Reynolds' 
equation  that  consumes  computer  time.  Solution  of  the  dynamics  equation 
determines  whether  the  response  to  a  disturbance  in  a  particular  degree 
of  freedom  will  grow  or  decay.  If  it  grows  the  system  is  unstable;  if  it 
decays  the  system  is  stable.  The  main  advantage  of  the  "Step-Jump"  method 
is  expeditious  determination  of  stability  thresholds.  Using  this  method 
(2)  a  three  shoe  bearing  was  analyzed  in  a  system  that  permitted  a  total 
of  5  degrees  of  freedom.  These  included  the  x  and  y  translations  of  the 
shaft  l'A  the  pitch  of  each  shoe.  These  are  the  predominant  motions  for 
any  tiiting-pad  journal  bearing  system. 

Examples  of  the  results  of  the  analyses  are  shown  on  Figure  11.  The 
bearing  system  is  stable  if  the  shoe  inertia  and  shaft  mass  are  such  that 
the  operating  point  is  below  the  stability  threshold  shown  on  Figure  11. 

This  plot  is  for  a  specific  compressibility  parameter  A,  length  to  diameter, 
L/D  ratio,  and  various  valves  of  pivot  film  thickness  Hp.  The  analysis  is 
confined  to  one  three-pad  bearing  with  fixed  pivots  operating  in  the  con¬ 
centric  position.  This  is  considered  a  representative  model  for  the  high 
speed,  low-load  conditions  encountered  by  miniature  cryogenic  turbomachinery. 
The  curve  indicates  that  as  the  shaft  mass  increases,  sc  does  the  stable 
range.  This  is  contrary  to  a  whirl  type  instability  i;.  which  an  increased 
mass  aggravates  the  whirl  problem.  Thus  the  indications  are  that  the  in¬ 
stability  is  due  to  pad  flutter. 

For  high  A,  variable  grid  geometry  was  again  utilized  to  avoid  numerical 
instability.  All  other  aspects  of  the  analysis  were  identical  to  the  con¬ 
stant  grid  situation. 

Appendix  III  presents  all  dynamic  performance  curves  such  as  that 
shown  on  Figure  11  for  the  range  of  variables  described  in  Section  1. 

Section  4  describes  procedures  for  determining  stability  of  a  particular 
bearing  configuration. 
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2.5  Limitations  of  the  Analysis 


There  are  a  number  of  restrictions  of  the  analysis  of  which  potential 
users  should  be  aware.  Removal  of  these  limitations  would  increase  com¬ 
plexity  and  make  the  production  of  design  information  impractical.  Most 
of  the  restrictions  are  concerned  with  secondary  effects  and  do  not  mate¬ 
rially  reduce  the  usefulness  of  the  design  information. 

a.  Pivot  friction  is  assumed  to  be  zero.  Naturally 
pivot  designs  should  be  as  frictionless  as  possible. 

b.  The  pivot  point  is  assumed  to  be  located  at  the 

•  surface  of  the  pad.  Thus,  no  account  is  taken  of 
the  surface  viscous  drag  moment  about  the  pivot 
which  is  usually  negligible.  It  is  good  practice 
to  locate  the  pivot  as  close  to  the  pad  surface  as 
practical.  (See  Section  4.) 

c.  The  dynamic  stability  thresholds  were  developed 
for  a  model  consisting  of  a  three-pad  bearing  in 
the  concentric  position  with  pad  motion  restricted 
to  pitch  mode  only.  This  is  the  primary  mode.  As 
discussed  in  Section  4,  the  thresholds  can  also  be 
applied  to  loaded  bearings  (non-concentric)  by  con¬ 
sidering  individual  pads  rather  than  a  complete 
bearing. 

d.  Unbalance  and  other  external  excitations  are  not 
considered  in  the  dynamic  analysis.  The  methods 
treat  self-excited  instabilities  due  to  fluid  film 
forces  only.  Normal  balance  limitations,  however, 
have  proven  more  than  adequate  for  gas-bearing 
miniature  cryogenic  turbomachinery,  although  operation 
should  not  be  in  the  vicinity  of  critical  speeds 
computed  on  the  basis  of  steady-state  bearing  stiffness. 
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3.  EXPERIMENTAL  PROGRAM 


A  decision  was  made  in  the  early  stages  of  the  project  to  utilize 
a  large  tilting-pad  journal  bearing  test  rig  (4-inch  diameter  shaft) 
that  was  available  at  FIRL  from  previous  programs.  The  reasons  for 
using  this  rig  rather  than  going  directly  to  the  .small  (1/2-inch  diameter) 
cryogenic  test  rig  used  on  the  previous  studies  (Reference  1)  were  as 
follows: 

a.  The  rig  has  a  prior  history  of  providing  excellent 
capabilities  for  conducting  experimental  tests  on 
tilting  pad  journal  bearings. 

b.  The  larger  size  permits  much  more  accurate  and 
easily  obtainable  measurements. 

c.  The  rig  incorporates  a  test  chamber  that  will 
allow  variations  in  the  environment  to  get  up  to 
the  high  lambda  values  necessary  for  steady  state 
testing. 

d.  It  could  be  made  available  for  testing  with  no  major 
modifications  or  major  preliminary  check-out  work. 

The  usual  operating  conditions  for  this  rig  are  shaft  speeds  up  to 
18,000  rpm. 

3.1  Description  of  Test  Rig 

Figure  12  shows  the  test  rig  mounted  on  its  wooden  base  for  atmo¬ 
spheric  testing.  Behind  the  rig  is  the  environmental  chantoer  that  was 
employed  for  high  A  tests.  The  rig  was  inserted  in  the  chamber,  which 
was  subsequently  sealed  and  evacuated  for  high  A  testing. 

The  rig  contains  two  bearing  pedestals  that  support  the  tilt-pads. 
The  pedestals  contain  provisions  for  mounting  of  capacitance  probes  for 
monitoring  the  shaft  and  selected  pads.  A  mylar  foil  supported  by  pulley 
mechanisms  is  used  to  apply  loads  to  the  shaft.  The  loading  mechanism 
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for  atmospheric  testing  employed  calibrated  dead  weights.  Note  that 
the  loading  mechanism  lifts  upward  on  the  shaft  and  is  designed  to 
remove  the  uead  weight  shaft  load  from  the  bearing.  The  total  shaft 
weight  is  92  lbs.,  so  the  load  range  on  each  bearing  can  vary  between 
0  and  46  lbs.  Figure  13  shows  the  test  rig  from  the  opposite  side 
and  gives  a  better  view  of  the  turbine  drive  and  thrust  bearing.  Also 
indicated  is  the  external  pressure  lines  that  feed  the  bearing  pads 
through  the  pivots  for  start-up  and  shut-down. 

Figure  14  shows  the  test  rig  modified  for  high  A  testing.  Since 
this  testing  is  accomplished  inside  the  environmental  chamber,  it  is 
cumbersome  to  use  dead  weight  loading.  A  bellows  load  cell  was  con¬ 
structed  and  loading  applied  through  air-pressure.  Calibrated  strain 
gages  were  used  to  establish  applied  loads.  Figure  15  shows  the  tank 
sealed  up  with  the  rig  inside,  ready  for  vacuum  environmental  (high  A) 
testing. 

3.2  Instrumentation 

Instrumentation  was  utilized  to  determine  the  following  measurements: 

a.  X  and  Y  translatory  motions  or  alternatively 
orbit  traces  of  the  shaft  at  each  journal 
bearing  location. 

b.  Relative  pad  pitch  to  shaft  motions  of  one 
fixed  pivot  pad. 

c.  Absolute  pitch  motions  of  one  fixed  pivot  pad. 

d.  Gross  shaft  displacement. 

e.  Shaft  speed. 

f.  Shaft  loads. 

g.  External  pressure  to  pad  pivots  (start-up, 

shut-down  and  operacion  below  10,500  rpm) . 

On  Figure  15  is  shown  the  strain  gage  read-out  bridge  used  for  de¬ 
termining  load  applied  by  the  pressurized  bellows.  The  pressure  gage  is 
fed  through  a  six-way  valve  and  is  used  for  establishing  pivot  pressure 
in  each  of  the  six  pads.  Mounted  above  the  tank  are  four  Wayne-Kerr 
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Figure  14.  Test  Rig  with  Pressurized  Bellows  Loading  Device 


Figure  15.  Test  Rig  Inside  Tank  for  High  A  Testing 
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distance  meter  boxes  that  amplify  capacitance  probe  signals  and  convert 
into  voltage  output  for  driving  dual  beam  oscilloscopes.  The  oscilloscopes 
measure  individual  vibrations  monitored  by  the  probes  or  the  probe  outputs 
can  be  combined  to  provide  orbit  traces.  Also  probe  outlets  are  fed  into 
a  Wayne-Kerr  electronic  micrometer,  which  measures  gross  displacements  of 
the  shaft.  The  electronic  micrometer  and  oscilloscopes  are  shown  on 
Figure  16. 

3.3  Pad  Configurations 

Figure  17  depicts  a  cross-sectional  drawing  of  the  two  fixed  pivoted 
pads.  A  photograph  of  the  two  pads  are  shown  on  the  bottom  of  Figure  17. 
The  inlet  orifice  and  surrounding  recess  are  clearly  shown.  The  variance 
in  pad  configuration  was  necessary  to  permit  installation  of  all  three 
pads  around  the  circumference  of  the  bearing.  Probes  for  measuring  abso¬ 
lute  pitch  motion  are  also  indicated.  Note  that  the  angular  distance 
between  pads  are  not  equal  for  the  test  rig.  This  differs  from  the  rec- 
on£j|nded  configuration  of  the  bearings  considered  for  miniature  cryogenic 
machines  and  for  which  design  information  has  been  generated.  However, 
the  objective  of  the  testing  was  to  validate  the  analytical  tools,  and  for 
this  purpose  a  special  analysis  of  the  test  bearing  was  completed. 

To  accomplish  high  A  stability  testing  it  was  necessary  to  produce 
sets  of  light  pads  since  it  was  anticipated  that  the  existing  heavy 
pads  would  go  unstable  prior  to  reaching  the  high  A  condition.  The  two 
different  pads  are  shown  on  Figure  18.  The  buttons  on  the  pads  are  used 
as  monitor  surfaces  for  the  capacitance  probes  for  measuring  absolute 
pitch,  roll  or  yaw.  When  the  light  shoes  were  installed,  relative  pitch 
motions  between  the  pad  and  shaft  were  measured  by  physically  attaching 
the  probes  to  the  pad. 


3.4  Difficulties  Encountered  During  Testing 

Two  major  problem  areas  occurred  during  testing.  The  first  was  con¬ 
cerned  with  the  tilting-pad  thrust  bearing  originally  installed.  This 
was  a  ■self'-equalining  Kingsbury  type  thrust  bearing  with  cylindrical  pivots 
used  on  the  leveling  plates.  The  pad  orientation,  when  operating  in  the 
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SHOE  CONFIGURATION  NO.  1 


SHOE  CONFIGURATION  NO.  2 


Figure  17.  Test  Bearing  Configuration 
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Figure  18.  Light  and  Heavy  Pads 
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horizontal  altitude  was  such  that  diverging  films  occurred  on  the  pad 
leading  edges.  This  caused  rubbing  contact  to  occur  between  the  shoes 
and  the  runner,  whicn  caused  friction- induced  whirl  and  increased  journal 
bearing  orbits  at  the  thrust  end.  To  alleviate  this  situation  a  new  com¬ 
pliant  surface  thrust  bearing  was  constructed.  The  bearing  consists  of 
separate  pads,  each  pad  supported  on  a  rubber  mount.  Deformations  of  the 
rubber  produce  some  equalizing  action  and  tilting  of  each  pad  to  form  the 
necessary  converging  wedge  for  hydrodynamic  load  capacity.  The  absence 
of  having  to  support  the  pads  on  pivots  eliminated  any  tendency  for  pad  hang¬ 
up  and  very  smooth  and  satisfactory  operation  resulted.  Figure  19  shows 
orbit  traces,  at  both  journals  with  the  original  thrust  bearing  installed. 
Figure  20  depicts  traces  with  the  compliant  surfaces  installed.  The 
smaller  orbits  and  absence  of  sub- synchronous  concepts  are  clearly  evident 
with  the  revised  thrust  bearing. 

The  other  major  difficulty  occurred  during  testing  to  determine  the 
effects  of  varying  unbalance.  Unbalance  response  results  are  shown  in 
Figure  21  for  a  0.0884  gram  unbalance  weight  located  at  the  turbine  end 
of  the  shaft  on  a  2-inch  radius.  At  18,000  rpm  (A  *  3.5)  the  unbalance 
force  is  5.59  lbs.  Comparison  of  scope  traces  with  those  for  the  balanced 
rotor  for  similar  operating  conditions,  shows  no  apparent  effects  of  the 
unbalance.  The  unbalance  weight  was  then  Incr .?  '.ed  to  0.4508  grams  (5X) 
with  no  significant-  effects  apparent.  (See  Figure  22.)  Also  shown  on 
this  figure  are  results  for  a  gross  unbalance  of  1.56  grams.  The  speed 
was  not  taken  past  13,500  rpm,  as  dangerously  high  responses  were  in 
evidence.  The  rotor  was  decelerated  for  shut-down.  At  6,500  rpm  excess 
shoe  flutter  caused  bearing  contact.  The  unit  did  not  seize  but  damage 
was  sufficient  to  warrant  re-surfacing  of  the  shaft  and  bearings. 

A  minor  problem  occurred  during  check-out  of  the  rig  with  the  light 
shoes  installed.  Attempts  to  balance  the  shaft  directly  in  the  rig  were 
not  successful.  It  was  necessary  to  use  a  conventional  machine  to  balance 
the  shaft.  Balancing  could  be  accomplished  with  the  heavier  shoes  in¬ 
stalled  (which  also  had  greater  pivot  friction).  The  significance  is  that 
the  lighter  shoes  appeared  to  respond  much  more  sensitively  to  unbalance 
than  the  h«-auier  shoas. 
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WEST  JOURNAL  BEARING  EAST  JOURNAL  BEARING 

TURBINE  END  THRUST  BEARING  END 


Figure  19.  Checkout  Runs  with  Pivoted  Thrust  Bearing. 

Sensitivies  are  lQOy-in/cm  and  5  m.sec/cm- 
Ful 1  Load. 
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Figure  20.  Checkout  Runs  with  Compliant-Surface  Thrust  Bearing 
Sensitivies  are  lOOp-in/cm  and  5  m.sec/cm-Full  Load 
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Shoe  Pitch  Motion 
Xi-D  a  Y2-A  Probe 
200  y-in./cm 
lm-sec/cm 


Shaft  Weight 


85  lb. 


Shaft  Orbit 
X5  &  Y5  Probes 
200  y-in./cm 
lm-sec/cm 
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Figure  21.  Unbalance  Response  of  Bearing  No.  1  (West  or  Turbine  End.) 

Unbalance  of  0.0884  Grams  at  Turbine  End  of  Shaft.  Speed 
is  18,000  rpm 
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Shoe  Pitch  Motion 
X5-B  &  Y4-C 
100  y-in./cm 
lm-sec/cm 


Shaft  Orbit 
X5  ft  Y4 
100  y-in./cm 
lm-sec/cm 


Speed  ■  18,000  RPM,  Unbalance  Weight  =  0.4508  gm. 


Unbalance  Weight  «  1.5608 


Figure  22.  Unbalance  Response  of  Bearing  No.  1  (West  or  Turbine 
End.)  Unbalance  Weights  of  0.4508  and  1.5608  Grams. 
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3.5  Steady-State  Test  Results 

Steady-state  verification  testing  was  accomplished  at  atmospheric 
conditions  (A  =  3.5)  with  the  heavier  set  of  pads  installed.  Attempts 
at  high  A  steady-state  testing  inside  the  evacuated  environmental  chamber 
with  the  hfeavy  pads  proved  difficult  and  dangerous  because  of  the  inherent 
tendency  to  produce  instabilities  as  loads  on  the  bearings  were  reduced. 

3.5.1  Comparison  Between  Theory  and  Experiment 

Figure  23  shows  the  comparison  between  theory  and  t'.st  for  a  single 
pad  steady-state  load  vs.  pad  film  thickness.  Test  results  on  the  full 
journal  were  converted  to  single  pad  results  for  convenience  of  comparing 
with  theory,  and  it  eliminates  another  variable,  namely  the  spring  pre¬ 
load.  As  indicated  on  Figure  23  comparative  results  agree  very  well  with 
theoretical  predictions  and  demonstrate  the  validity  of  the  governing  equa¬ 
tions  and  methods  of  solution. 

3.6  Dynamic  Test  Results 

The  most  significant  dynamic  testing  was  accomplished  with  the  light 
shoes  at  high  A  conditions ,  since  this  provides  the  best  environment  for 
pad  instability.  Dynamic  testing  with  the  heavy  shoes  at  low  A  (3.5)  was 
inconclusive. 

Table  1  is  a  summary  of  the  pertinent  dynamic  runs  conducted  with  the 
light  shoes.  Variables  included  the  bearing  pre-load,  ambient  pressure 
and  external  load.  Combined  pre-load  and  external  load  was  converted  to 
that  on  a  single  fixed  pad  to  facilitate  comparisons  between  theory  and 
experiment. 

Figure  24  shows  oscilloscope  orbit  traces  of  the  two  bearings  at 
18,000  rpm  at  atmospheric  pressure  (A  =  3.5).  The  orbit  size  for  the 
thrust  end  journal  is  approximately  twice  that  of  the  turbine  end  journal. 
The  conical  orbit,  which  was  characteristic  for  the  light  pad  installa¬ 
tion,  could  be  due  to  slight  residual  unbalance  or  possibly  influence  of 
the  thrust  bearing. 
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Figure  23.  Comparison  of  Theory  and  Test  for  Steady-State  Pad  Load 
vs.  Pivot  Film  Thickness. 
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TABLE  1 


Summary  of  Results  of  Dynamic  Testing  of  Light  Pad  Bearings 


Shaft  speed  **  18,000  rpm 

M  -  2.61  X  10"9  lb. -sec. /in. 2 
Shaft  diameter  *  4  inches 
C  =  0.0015  inch 


Pad  size  =  94.5  degrees,  3.3  inches  long 
Pad  weight  =  2.2  lbs. 

Shaft  weight  *  92.5  lbs. 


Pad 

Inertia 

-  0.0135  lb 

.-sec. 

2-in. 

Ambient 

Foil 

Pad 

Test 

Pre- 

Pressure 

Load 

Load 

No. 

Load 

psia 

lbs. 

lbs. 

1 

22  .r2 

15-0 

0 

53.2 

2 

20 

42.2 

3 

30 

36.7 

4 

40 

31.2 

5 

45 

28.5 

6 

12.0 

0 

53.2 

7 

10 

47.7 

8 

20 

42.2 

9 

30 

36.7 

10 

40 

31.2 

11 

8.75 

0 

53.2 

12 

=10 

47.7 

13 

7.5 

0 

53.2 

14 

6.25 

0 

53.2 

A  Remarks 
3.5 


Threshold  between  45  -  50 
lbs.  foil  load.  Pad  load 
=  27  lbs. 

4.37 


Threshold  between  40  -  45 
lbs.  foil  load.  Pad  load 
“  30  lbs. 

6.0 

Threshold 

7.0 

8.4  Threshold  at  P  =6  psia 

a 

@  53.2  pad  load,  rub  occured 
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Test 

.Bis. 

Ujgd 

fcMent 

Pressure 

D*fa 

Foil 

Load 

lbs. 

Fad 

Loa4 

tbs. 

A 

Remarks 

15 

36.5 

15  A 

0 

64.4 

3.5 

16 

10 

58.9 

17 

25 

50.6 

IS 

40 

42.4 

19 

11.75 

0 

64.4 

4,47 

20 

8.75 

6,0 

21 

7.5 

7.0 

Threshold 

22 

7.9 

15,0 

0 

42.  r 

3.5 

23 

10 

36.6 

24 

13 

->3.9 

25 

20 

31.1 

26 

12.5 

0 

42.1 

4.2 

27 

10.0 

5JJ5 

Threshold  @  ambient 

9.5  psia. 
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West  Shaft  Orbit 
Turbine  End 

Sensitivity  is  300  y  in. /cm. 
Orbit  Size  is  «  100  y  in. 


East  Shaft  Orbit 
Thrust  Bearing  End 
Sensitivity  is  100  y  in. /cm. 
Orbit  Size  is  =  200  y  in. 


Shaft  Orbits  at  18,000  RPM; 

Ambient  Pressure  is  15  psia; 

Bearing  Pre-Load  is  22.2  lbs.,  (c7c  =  .92); 

Pad  Inertia  is  .0135  lb. -in. -sec. 2 

Figure  24.  Bearing  Shaft  QrLit,  18,000  rpm,  15  psia,  22.2  lbs 
Preload 
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Figure  25  shows  the  effect  of  varying  bearing  load  and  pre-load  on 
the  thrust  end  bearing  (higher  orbit  end) .  Vertical  columns  have  constant 
pre-load  in  descending  order  left  to  right.  Bearing  load  increases  up¬ 
ward.  In  the  first  column  a  relatively  heavy  pre-load  is  applied  (36.5 
lbs.).  The  bearing  showed  no  signs  of  distress  until  the  bearing  load 
was  reduced  to  some  9  lbs.  (not  shown)  when  evidence  of  instability  was 
detected.  When  the  pre-load  was  reduced  to  22.2  lbs.,  instability  set  in 
at  a  bearing  load  of  approximately  13  lbs.  A  further  reduction  in  pre-load 
to  7.9  lbs.,  resulted  in  an  instability  setting  in  at  approximately  30  lbs. 
load. 

Figure  26  shows  the  effects  of  variations  in  ambient  pressure  and 
pre-load.  Bearing  load  was  maintained  at  a  constant  value.  Lowering  the 
ambient  pressure  reduces  the  stable  regime.  It  was  not  possible  to  obtain 
photographs  of  the  instabilities  for  each  situation  because  it  occurs  very 
suddenly.  Increasing  pre-load  to  a  limit  is  beneficial,  but  if  pre-load 
is  excessive  (overloaded  bearing)  instability  again  sets  in. 

Thus,  summarizing  it  can  be  said  that: 

1.  Instabilities  set  in  very  rapidly. 

2.  Increasing  pre-load  is  usually  beneficial. 

3.  Decreasing  ambient  pressure  (increasing  A) 
reduces  stability. 

3.6.1  Comparison  Between  Theory  and  Experiment 

Figure  27  depicts  stability  threshold  as  a  function  of  pad  load  coef¬ 
ficient  and  compressibility  parameter  A.  Experimental  points  fall  reason¬ 
ably  close  to  predicted  values  except  at  the  high  A  condition,  where  the 
theoretical  curve  is  conservative. 

Reasons  for  variations  between  experiment  and  theory  are  as  fellows: 

a.  The  analysis  does  not  account  for  unbalance 
response.  It  only  examines  self-excited 
instabilities. 
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Bearing  Load  =  39.2  lbs. 


Bearinc  Load  =  32  lbs. 


Bearing  Load  «  39.2  lbs. 


Ambient  Pressure  =  15  psia 


Ambient  Pressure  «  15  psia 


Ambient  Pressure  »  11.8  psia 


aT7\  7 


P,  =  12  psia 

a 


Ambient  Pressure  *  12.5  psia 
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b.  The  dynamic  analysis  is  a  small  perturbation 
theory  and  cannot  predict  bounded  orbits. 

c.  The  effects  of  pivot  friction  are  not  accounted 
for. 

d.  In  the  test  rig,  pivot  film  thicknesses  are  not  equal  on  all 
pads  at  zero  applied  load,  as  is  assumed  in  the  analysis. 
This  is  significant  at  high  A  where  the  analysis 

predicts  that  stability  is  sensitive  to  pivot 
film  thickness.  (See  Appendix  III). 
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4.  DESIGN  PROCEDURES 


The  intent  of  this  section  is  to  present  general  computational 
methods  for  establishing  geometrical  and  performance  parameters  using 
the  design  information  included  in  the  appendices.  A  complete  and 
all-inclusive  design  guide  was  beyond  the  scope  of  the  program,  but 
this  report  contains  adequate  design  information  and  sufficient  pro¬ 
cedural  discussion  to  enable  designers,  reasonably  well  versed  in  gas 
bearing  technology,  to  produce  a  successful  design.  In  addition,  a 
practical  design  configuration  which  has  had  good  success  at  FIRL  is 
suggested.  The  computations  are  divided  into  two  major  areas: 

(a)  establishing  steady-state  performance  and  (b)  determining  whether 
the  selected  system  is  stable.  Much  of  the  material  presented  below 
has  been  extracted  from  Reference  1  with  some  minor  modifications  and 
extensions  for  off-design  operation. 

4.1  General  Procedure  for  Establishing  Steady-State  Performance 

1.  Given:  For  any  problem  the  following  information  is  generally 

given. 

a.  ambient  pressure,  p  ,  psia 

b.  shaft  speed,  N,  rpm 

c.  shaft  diameter,  D,  in. 

d.  maximum  bearing  load,  W,  lbs. 

2 

e.  lubricant  viscosity,  p,  lb. -sec. /in. 

2.  Primary  considerations  for  determining  the 

length  of  the  bearing  are  maximum  load,  startup 

and  roll  stability.  At  startup  the  bearing 

is  only  subjected  to  gravity  load. 


W 

2RL 


<_  3  psi 


for  startup  where  W 


start-up  load  on  bearing 
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■r~T  <  10  psl  for  max  loading  when  operating 
/RL 

A  general  rule  of  thumb  is  to  limit  L/D  <  1.5 
to  prevent  roll  dynamic  problems. 

3.  Initially  investigate  performance  at  a  compress¬ 
ibility  parameter  A  =  1.5,  in  order  to  optimize 
the  design  for  minimum  power  loss  (3) .  Establish 
machined  in  clearance  C  from: 


4.  Compute  load  coefficient  for  three  shoe  bearing. 


„  W 

C  “ - 

LT  p  RL 
*a 


5.  Enter  the  appropriate  full  bearing  design  chart  in 

Appendix  I  at  the  computed  value  of  C  .  For  vary- 

LT 

ing  values  of  the  pre-load  factor  C'/C  tabulate  the 
following : 

a.  Bearing  eccentricity  ratio,  t'. 


will  always  equal  for  the  geo  ’"tries  con¬ 
sidered.  When  shaft  displacement  is  between 
pads,  Hpl  and  will  be  less  than  When 

displacement  is  directly  into  pad  3,  will 
be  less  than  Hp^  and  Hp2. 


c.  Proper  interpretation  of  the  curve  requires 
the  following  procedure: 

1.  The  intersection  of  the  value  of  C  with 

LX 

a  selected  C'/C  curve  establishes  e'. 


At  the  established  t'  go  to  the  Hp^  2  and 
Hp3  curves  at  the  corresponding  values  of 

C'/C  and  read  off  the  pivot  film  thicknesses 
from  the  right  hand  ordinate.  Compute  di¬ 
mensional  pivot  film  thicknesses  from: 


=  H  ,  X  C 
Pi 


i  -  1,  2,  3. 
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6.  Determine  the  non-dimensional  minimum  film  thickness 
for  each  pad  from  the  appropriate  minimun  film  VS. 
pivot  film  thickness  curve.  Then  dimensionalize  by: 

h  =  H  X  C. 
m  tn 

7.  Determine  the  pad  friction  coefficient  F  for  each 

pad  from  the  appropriate  F  vs.  curve.  The  friction 

force  for  each  pad  F^  is  determined  from: 


F,.  -  FP  CL 
fi  a 


i  =  1,  2,  3 


and  the  friction  power  loss  in  watts  is  determined 
from: 

Tl?  V  D 

»f  ■  w  * x  746- 

8.  From  appropriate  individual  pad  load  coefficients 
vs.  pivot  film  thickness  curve  establish  pre-load 
on  spring-loaded  shoe: 

W  =  C  F  RL.  (VJ  =  load  or  preloaded  pad) 

Li  cl 

9.  From  the  above  procedure  major  operating  performance 
at  the  design  condition  is  established.  With  regards 
to  selection  of  the  pre-load  setting  C'/C,  there  are 
a  number  of  constraints  that  are  applicable. 

a.  C'/C  is  the  ratio  of  the  pivot  circle  clearance 

to  the  machined  clearance  of  a  pad.  A  small  value 
implies  tight  film  thickness  with  a  consequently 
large  power  loss.  A  large  value  may  unload  the 
shoe  excessively  with  possibilities  of  shoe  vi¬ 
bration.  Also,  large  values  of  C'/C  can  also  cause 
a  phenomenon  known  as  leading  edge  lockup,  where 
the  leading  edge  of  the  shoe  is  forced  into  the 
rotating  shaft.  As  a  general  rule: 

0.3  £  C'/C  _<  0.8. 

In  no  instance  should  C'/C  exceed  unity. 

b.  If  the  bearing  design  is  tc  incorporate  a  spring- 
loaded  shoe,  then  the  pre-load  factor  C'/C  is 
directly  related  to  the  spring  pre-load.  For 

a  spring-loaded  shoe  C'/C  equals  the  pivot  film 
thickness  H^.  From  this  the  load  on  the  spring 
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leaded  shoe  can  be  determined.  The  spring 
design  should  be  such  that  normal  shoe  dis¬ 
placements  will  cause  negligible  change  in 
the  spring  pre-load,  or  C'/C,  and  the  bearing 
will  behave  as  predicted. 

10.  Other  constraints  are  power  loss  and  minimum  film 
thickness.  Maximum  power  loss  is  dependent  upon  the 
application,  and  allowable  minimum  film  thickness 
depends  upon  the  size,  production  accuracy,  and  di¬ 
mensional  changes  caused  by  structural  and  thermal 
deformations.  For  shaft  sizes  1/P.  inch  and  less  a 
suggested  absolute  minimum  is  0.00025  in.  and  for 
shafts  up  to  1  inch  a  suggested  minimum  is  0.0005  in. 

11.  Once  a  design  is  selected,  it  is  often  desirable 
to  determine  off-design  operation.  The  procedure 
varies  depending  upon  whether  the  pad  ii:  spring 
loaded  or  whether  the  pivot  circle  is  considered 
fixed  (three  fixed  pads). 

a.  Spring-loaded  shoe  -  Over  the  shaft  displace¬ 
ment  range  the  spring  pre-load  variation  is 
insignificant,  and  thus  the  load  on  the  spring- 
loaded  shoe  is  constant.  A  constant  load  im¬ 
plies  a  constant  pivot  film  thickness  H^,  and 

off-design  performance  is  obtained  by  following 
lines  of  constant  on  the  appropriate  full 

bearing  load  vs.  eccentricity  plots.  The  curves 
in  Appendix  I  are  segregated  into  fixed  pivots 
on  all  three  pads  and  one  spring-loat  <sd  pad  with 
two  fixed  pivoted  pads.  For  most  cryogenic  ap¬ 
plications  a  spring-loaded  pad  is  recommended 
because  of  strong  thermal  and  centrifugal  growth. 

b.  Fixed  pads  -  If  all  pads  are  fixed,  and  thermal 
or  centrifugal  shaft  growth  will  not  change  the 
original  pivot  circle  relationship,  then  the 
pre-load  varies  and  performance  is  obtained  by 
following  lines  of  constant  C'/C  on  the  proper 
fixed  pivot  performance  curves. 

4.2  Sample  Problem  -  Steady-State  Performance 

As  a  sample  problem  the  design  of  a  typical  bearing  for  a  miniature 
cryogenic  turbomachine  will  be  described.  The  following  information  is 
given : 
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1.  Shaft  Diameter,  D  =  0.5  inches 

2.  Ambient  pressure,  pa  =  15  psia 

3.  Shaft  speed,  N  -  200,000  RPM 

4.  Fluid  temperature,  T  =  100 °K 

5.  Bearing  radial  load,  W  =2.5  lbs.  (maximum,  primarily  due  to 
unbalance  loading) 

-9 

6.  Lubricant  is  gaseous  nitrogen  with  viscosity,  jj  =  1.01  X  10 

2 

lb. -sec. /in. 

a.  bearing  length 


W 

2RL 


<  10 


L 


min 


W 

20R 


2.5 

20  <|) 


0.5  inches. 


This  length  was  used  as  the  first  approximation,  and 
was  found  to  be  satisfactory  in  all  respects. 


b.  machined  in  clearance 


C  =  7.27  X  10  4  A~1/2 


To  optimize  the  bearing  for  viscous  power  loss,  A  =■  1.5. 
Dimensional  inaccuracies,  thermal  growth  of  the  shaft, 
etc.  would  result  in  reduction  of  the  tip  pad  clearances, 
C;  therefore,  a  A  =  3.5  was  also  considered  to  provide 
performance  estimates  over  the  possible  clearance  range. 

A  C 

1.5  5.937  X  10~4  in. 

3.5  3.887  X  10~4  in. 

c.  load  coefficient 


C 


LT 


W  _  2.5 _ 

P  RL  15  X  0.25  X  0.5 
a 


1.333 


d.  following  the  procedures  of  steps  5,  6,  7  of  paragraph 
4.1.  Table  2  can  be  produced.  The  numbered  columns 
are  generated  as  described  below: 


49 


mUMd  . . . . . ,  1 . rYiiifMiiiiiiWiii 


o 

O 

u> 

li 

CQ 

•  A 

o 

O 

o 

II 

0 

*  A 

<M  . 

c 


u 

0) 

t/» 

I 

-Q 


ON 

I 

O 


O 

t  ^ 

Clt  r— 
Q.  O 

■ 

cn  r— 


CM 


S*. 
<t5 
QJ 
LlJ  CD 
_J 

CD  r- 
<C  05 


O 

T? 

fO 

Q. 

o> 

c 


fO 


CO 

O. 


LO 


II 


LO 

CM 

II 

+J 

3: 


in 

at 

o 


(O  tn 

CL 

ro 

•  « 

s:  ii 
a. 

cc  c 
CD  *« 

o  in 
o  in 


o 

o 

CM 


s 

X 

-©■ 


m 

o 

n 

Q 


in 

o 

ii 


^  * 
<2  o  c 

(12) 

0.59 

0.77 

0.80 

0.83 

m  CO  rH  vo 
rt  -3  UN  UN 
•  •  •  . 

<t  • 

C0  O  c 

a  rH  -rl 

.c  x 

(11) 

H  ON  in  H 

rx  oo  on  o  i 

•  •  •  •  | 

O  o  O  *H  1 

IO  rH  00  CM 

-■t  tn  tn  vo 

•  •  •  • 

CM  si- 

•»  o  • 

H  H  13 
O.  X  -H 

rC 

(10) 

2.33 

3.15 

3.96 

5.  71 

co  co  eg  o 
tn  o  v©  eg 
•  •  •  • 
eg  co  oo  <f 

Pf 

watts 

/■S 

o> 

w 

CM  «3-  m  s3-  1 
•  .  .  •  j 

o  oo  rx  vo  j 

•4*  co  »h  in 

•  .  «  . 

rt  c|  h  O 
H  r;  H  H 

CO  5 

sTS 

«  a 

•H 

-c-N 

OO 

'w' 

•j  voon 
co  m  cm  ix  i 

•  •  .  •  j 

S3"  CO  CO  CM  | 

cm  tn  ro  oo 

IX  CM  IX  ^ 

•  •  •  * 
IO  ITI  sf  sr 

CO 
O 
U-4  rH 

r\ 

r^. 

w 

4.34 

3.56 

3.20 

2.73 

eg  in  co  oo 

N  OJ  N 
•  •  •  • 
tn  m  <j-  <f 

co 
Pn  ■ 


VO 


co 


CM 

tH 

.ca 


co 


<f 

I 

o 


co 

<r» 

• 

m 

i 

u 

•  « 
4J 

£ 

•rl 

a< 

o 

u 

e 


T3 

a) 

4-1 

o 

at 

u 


in 


ix  m 

<3N  CM  00  Mf  I 
•  •  •  •  j 
co  oo  CM  CM  | 


m 

CM  ON  IX 

ix  \o  in  m  i 
•  •  .  •  I 

rl  rt  H  rl  j 


CM 

On  MT  CO 

H  rx  VO  <r 


CM 

cm  m  no  ix  rx 
H  rl  H  r)  H 


CM  ix  CM 
ON  CO  VO  vO  CM 
co  m  vO  O'  H 


m  in  o  in 
rx  cm  co  no  on 
m  vo  vo  ix  rx 


co  sr  m  rx  oo 


sr 

l 

o 


lx 

oo 

00 

• 

CO 

II 

(J 


4-1 

o 

> 

•H 

Pm 

O 

4-1 

0 

•H 

at 

4J 

o 

0) 

M 


T3 

tfl 

O 

rJ 

in 

• 

co 

it 

< 


m  o  o  m 

00  CM  in  rH 
•  is3  •  • 

rx  r^  vo  vo 


m  o  O  m 

O'  IX  vj-  CM 
•  •  •  • 
Mf  <J-  If  <t 


tn  in  in  m 

xr  CM  o  ON 
H  HH  HO 


"}SOo 
rl  Cl  in  vo 


m  oo  co  oo 
NO  lx  on  o 


rx  m  o  cm 
vo  h  in  oo 
m  vo  vo  vo 


m  vo  rx  co 


-3- 

o 


lx 

CO 

ON 

• 

m 


to 

4J 

£ 

•H 

&4 

e 

4) 

01 

3 

0) 

pq 

•O 

tfl 

>3 

m 


«3-  co  m 
M3  oo  O' 


XO  O  -O 
h  in  ix 
•  •  . 
CM  Cl  VO 


H  CM  O 
Ix  ON  H 


OO  sr  ON 

•  •  • 
HOB 


CM  rH  O 
OM00 
•  •  * 
in  -o  co 


CM  rH  3 
O  Mf  00 
•  «  » 
m  -a-  co 


H  no  H 
UN  ON  -J 
•  •  • 
Mt  CO  CO 


rx  cm  ix 
O  IX  CO 


rH  O  O 


CM  CM  CM 
rx  vo  m 


in 

NO  ON  co 
co  m  rH 


m  m 
cm  in  oo 

H  H  H 


O  CO 
CM  vO  ix 
00  O'  CM 


cm  co  m 


rx 

oo 

oo 

• 

CO 

II 

O 


> 

tH 

04 

c 

a) 

ai 

3 

at 

(O 


o 

»3 

m 

co 

n 

*3 


BOO 
vo  oo  oo 

<f  vO  ix 


m  co 
OnOnO 
CM  C3N  ON 
•  •  • 
tH  H  CO 


<f  fs. 

vo  h 
in  !>•  on 


o  o  o 
•  •  • 

O  vo  CO 
NH  H 


O'  H  ^ 

<rcom 
#  •  • 
oo  vo  tn 


5n  H  s} 

<t  oo  m 


oo  vo  in 


m  m 
vo  co  vO 
•  •  • 
H  On  h» 


in  tn 
vDQO  O 
•  *  • 
eg  h  H 


o  m  o 

m  n  n 


>4-  co  co 


m 

H  o  eg 
co  in  o 


m  m 
<5*  o  co 
rH  eg  eg 


vocosr 
m  >o  o 


eg  co  in 


50 


(1)  From  appropriate  Figure  of  vs. 

e'  at  intersection  of  and  C'/C 
curve . 


(2), (3)  Maintain  e'  and  read  from 
intersection  of  z'  with  C'/C 


(4),  (5)  From  figure  of  friction  power  loss 

(6)  EF  =  Fx  +  F2  +  F3  =  2FX  +  F3 

(7)  F  =IFXP  CL 

f  a 

(8)  M  =  Ff  X  R 


(9)  pf  = 

(10)  , (11) 


Mf  X  N 
63,000 


X  746 


h  =  H  X  C 
P  P 


(12)  h  =  H  x  C  where  H  from  appropriate  curve  of  H  vs. 
v  /  m  m  m  m 


H 

P 


The  minimum  film  thickness  and  power  loss  Vs.  pre-load  parameter  are 
plotted  on  Figures  28  and  29.  Note  that  the  minimum  film  thickness  approaches 
a  limiting  value  for  most  cases  considered,  and  an  optimum  value  of  C'/C 
is  about  0.5,  since  it  provides  nearly  maximum  film  thickness  and  reasonable 
powvjr  loss. 


4.3  Sample  Problem  -  Dynamics 

Given,  information  of  the  type  presented  in  Appendix  III,  the  de¬ 
signer  is  now  in  a  position  to  design  both  tilting  pad  bearings  and 
rotors  which  will  be  stable  over  the  desired  operating  range.  All  the 
important  parameters  such  as  speed,  ambient  pressure,  viscosity,  and 
shoe  and  shaft  mass  and  inertia  properties ,  as  well  as  bearing  dimen¬ 
sions  such  as  bearing  diameter,  length,  clearance,  pivot  location  and 
pad  arc,  can  be  specified  by  the  designer  to  assure  a  stable  system.  As 
an  example,  given  a  desired  operating  speed,  ambient  pressure  and  gas, 
the  designer  can  calculate  the  required  bearing  dimensions  and  mass, 
and  rotor  properties  that  will  result  in  a  stable  bearing-rotor  system. 

On  the  other  hand,  if  the  designer  has  the  bearing  and  rotor  properties, 
he  can  establish  the  threshold  of  instability  and,  consequently,  specify 
the  upper  speed  limit  above  which  the  bearing-rotor  system  is  unsatisfactory. 


51 


(SilVM)  *d  'SSOl  d3M0d  NOIXOIdd 


UJ 

or 

cl 


53 


Figure  29.  Friction  Power  Loss  vs.  Pre-Load  Parameter  at  Design  Conditions,  Sample  Problem 


As  an  example,  consider  Che  typical  cryogenic  tilting-pad  bearing-shaft 
combination  design  for  which  the  steady-state  analysis  has  been  described. 
The  pitch  moment  of  inertia  for  a  particulcr  pad  configuration  was 
given  as: 

1^-  .0205  X  10“4  y  lb.-sec.2-in. 

where 

y  »  0.283  lb. /in. 3  for  steel 
3 

*  0.162  lb. /in.  for  titanium. 

The  shaft  mass  is: 

M  *  1.2081  X  10“3  y  lb.-sec.2/in.4 
s 

Operating  conditions  for  which  it  is  required  to  determine  stability  are 
specified  as: 

N  *  200,000  rpm 
p  *  Ji  psia 

C  -  0.0003887  in. 

R  -  0.25  in. 

L  *»  0.50  in. 

To  enter  the  stability  curves  of  Figures  30  and  31,  it  is  necessary  to 

compute  the  dimensionless  mass  parameter  for  the  shaft,  M  ,  and  the 

s 

dimensionless  shoe  inertia  parameter,  1^.  These  computations  follow: 

M  CO2  (1.2081  X  10"3y) (.0003887)  (~  200,000)2 
\  *  4p  RL  “  (4)  (15)  (0.25)  (0.5)  3  0,2746  X  10  y 

Inserting  values  for  Y  one  obtains: 

M  *  7.77  for  a  steel  shaft 
__s 

M  *  4.45  for  a  titanium  shaft 
s 

and 

IDC  2  (.0205  X  10“N) (.0003887) 2000, 000) 2 

I  *  — — —  «  - .  738  Y 

P  4paR3L  (4) (15) (.25)3(.5) 
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3 

For  steel  shoes,  y  =  .283  lb. /in.  and 

I  *  .209 
P 

For  titanium  shoes,  I  =  .120 

P 

Figures  30  and  31  show  the  various  combinations  of  shoe  and  shaft 
materials  for  operation  at  A  =  1.5  and  3.5  respectively.  The  combina¬ 
tion  of  the  titanium  pads  and  steel  shaft  provide  stable  performance  at 
both  values  of  A  and  should  be  the  selected  combination. 

This  sample  problem  describes  one  possible  use  of  the  dynamic 
stability  design  charts.  Instead  of  material  combinations,  it  would 
have  been  also  possible  to  vary  the  inertia  properties  of  the  pads  or 
the  mass  of  the  shaft. 


The  limitation  of  concentric  operation  may  be  circumvented  by 
examining  what  occurs  on  a  single  pad  rather  than  the  entire  bearing. 

For  example,  for  a  heavily  loaded  situation,  the  lightly  loaded  pad 
would  have  a  high  pivot-film  thickness  H^,  and  its  dynamic  stability 
can  be  established  by  the  corresponding  threshold  curve.  Similarly 
for  the  loaded  pads.  On  Figure  30  the  titanium  pads  steel  shaft  com¬ 
bination  is  stable  for  both  high  and  low  values  of  pivot  film  thickness. 


11^,  so  that  acceptable  operation  o\er  a  wide  load  range  could  be 
anticipated. 


4.4  Practical  Design  Considerations 

A  useful  design  concept,  which  is  a  ramification  and  improvement  of 
a  design  which  has  been  successfully  employed  at  FIRL  is  shown  on  Figure 
32.  It  consists  of  a  split  housing  to  permit  locking  of  the  adjustable 
fixed  pivot  pins.  The  pivot  pins  themselves  are  accessible  through  ports 
in  the  outer  housing  sleeve.  The  spring  mounted  pivot  is  loaded  through 
a  beam  type  spring  which  can  be  pre-loaded  a  desired  amount  by  adjusting 
screws  and  pre-load  mount.  The  pivots  are  the  ball  and  socket  type  with 
the  socket  diameter  slightly  larger  than  the  pivot  diameter.  The  pivots 
are  located  close  to  the  face  of  the  pad  to:  (1)  keep  pitch  inertia  low; 
(2)  keep  lock-up  friction  moment  low;  and  (3)  minimize  camming  during 
initial  alignment. 
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Figure  30.  Dynamic  Stability  Curve  for  A*1.5,  Sample  Problem 
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Figure  31.  Dynamic  Stability  Curve  for  A=3.5,  Sample  Problem 


i 

\ 


Figure  32.  Journal  Bearing  Assembly 


In  general  it  is  desirable  to  include  the  following  capabilities  in 
design: 

a.  Have  a  cartridge  type  assembly  for  ease  of 
replacement . 

b.  Make  pivots  adjustable  for  flexibility  in 
pivot  clearance  setting  and  to  facilitate 
alignment  of  the  thrust  bearing. 

c.  Insure  that  pre-load  of  spring-mounted  pivot 
is  externally  adjustable. 

d.  Keep  pre-load  spring  soft  so  that  pre-load 

•  remains  essentially  constant  and  centrifugal 
and  thermal  growths  can  be  readily  accomodated. 

Successful  materials  include  tungsten  carbide  coated  surfaces 
applied  to  410  stainless  steel  base  material  against  stainless  or 
titanium  shafts.  It  is  also  possible  to  coat  the  pads  with  aluminum 
oxide  and  mate  with  a  journal  coated  with  nickel  chromium  bonded 
chromium-carbide.  A  flame  plated  tungsten  carbide  journal  against 
P692  carbon  graphite  pads  is  another  possibility. 
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6.  NOMENCLATURE 


=  machined  clearance,  in. 

*  pivot  circle  clearance,  in. 

=  preloaded  factor 

=  pad  load  coefficient  -  W/P  RL 

a 

=  bearing  load  coefficient  =  VJ/P  RL 
•  t  a 

=  diameter  (shaft),  in. 

=  dimensionless  time  step 

=  Shaft  displacement,  in. 

=  dimensionless  friction  =  F-/p  CL 

f  ra 

=  viscous  friction  force,  els. 

=  viscous  friction  force  on  shaft,  lbs. 

=  dimensionless  film  thickness  =  — 

c 

=  dimensionless  minimum  film  thickness  **  h  /C 

m 

=  dimensionless  pivot  thickness  =  hp/C 

*  film  thickness,  in. 

=  minimum  film  thickness,  in. 

=  pivot  film  thickness,  in. 

2 

=  pitch  moment  of  inertia  of  shoe,  lb-sec  -in. 

=  dimensionless  pitch  moment  of  inertia  of  shoe  = 
=  bearing  length,  in. 

=  dimensional  viscous  friction  moment,  in-lbs 


I  C  Sr 


4p  RL 
a 


=  shaft  mass,  lb-sec  /in. 

=  dimensionless  shaft  mass  = 
=  shaft  speed,  rpm 
=  normal  unit  vector 


M  C  S2 

J 

4p  RL 
ra 


*  pressure,  psia 

=  ambient  pressure,  psia 

15  total  bearing  power  loss,  watts 
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«  shaft  radius,  in. 

•  length 

•  temperature,  °K 

■  distance  in  length  direction 

-  time,  seconds 

*  velocity  vector 

-  volume,  in. 

«  bearing  or  pad  load,  lbs 

*  pad  angle,  rad 

*■  angle  from  load  vector  to  pivot,  rad 

3 

-  specific  weight,  lb/ in. 


*  bearing  eccentricity  ration  *  e’/C 

.  .  _  lb-sec2 

■  mass  density  - — 


«  circumferential  distance,  radians 
»  compressibility  parameter  * 

*  lubricant  viscosity,  lb-sec/in, 2 


-  angle  from  leading  edge  of  shoe  to  pivot 
•  shaft  angular  speed  rad/ sec 
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APPENDIX  I 

STEADY-STATE  ANALYSIS  AND  DESIGN  CURVES  FOR  TILTING  PAD  JOURNAL  BEARING 


1.  INTRODUCTION 

The  design  charts  for  pivoted  pad  bearings  (three  pads)  are  included 
in  this  appendix.  The  charts  have  been  developed  for  the  following 
parameters : 


Comp  res  s ib i li ty 
Number  A 

R/L 

1.5 

3.5 

0.25 

5.0 

0.50 

10.0 

1.0 

20.0 

The  pivot  location  <J>/a  =  0.65,  and  the  angular  extent  of  the  pad 
a  =  100°,  and  the  angular  distance  between  pivots  26  =  120°  are  constant 
for  all  information  shown. 

For  each  combination  the  following  design  charts  are  included: 

a)  Bearing  Load  Coefficient,  C  „  and  Pivot  Film  Thickness, 

Hp  vs.  Bearing  Eccentricity  Ratio  e'  with  or  without  a 
spring  loaded  pad  for  load  between  pads 

b)  Same  as  a)  for  load  into  pad 

c)  Friction  Coefficient,  F  vs.  Pivot  Film  Thickness,  H  , 

for  a  single  pad  ^ 

d)  Lift  Coefficient,  G^,  vs.  Pivot  Film  Thickness,  H  , 

for  a  single  pad  ** 
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e)  Minim  un  Film  Thickness,  vs.  Pivot  Film  Thickness, 

Hp  for  a  single  pad 

2.  PROCEDURE  FOR  GENERATING  STEADY-STATE  BEARING  PERFORMANCE  CURVES 

The  usual  procedure  for  determining  steady-state  performance  of  a 
gas-lubricated  tilt,  ng-pad  journal  bearing  is  to  develop  field  maps  for 
individual  pads  with  the  aid  of  a  steady-state  computer  program.  The 
pivot  position  is  output,  since  it  is  the  center  of  pressure  and  carmot 
be  known  until  the  pressure  distribution  is  formulated.  From  these  field 
maps,  three  pad  bearing  performance  can  be  obtained. 

For  every  combination  of  A  and  R/L,  the  generation  of  field  plots 
require  many  data  points  and  a  considerable  amount  of  computer  time.  In 
addition,  since  the  optimum  pivot  location  <p/ ct  is  0.65,  many  of  the  field 
plot  data  points  are  extraneous,  and  unnecessary.  Thus  a  new  scheme  was 
devised  that  was  much  more  efficient  in  the  production  of  pertinent 
steady-state  data.  Reference  (4)  describes  an  analysis  for  obtaining 
dynamic  performance  of  a  Tilting-Pad  Bearing  System.  The  method  dis¬ 
cretizes  time  into  small  finite  elements.  The  process  starts  from  a 
fixed  shaft  position  and  pivot  position  and  an  arbitrary  pitch  angle. 

Also  it  starts  with  an  assumed  pressure  distribution  on  the  pad.  The 
sequence  for  obtaining  steady-state  data  using  the  dynamics  program  is 
as  follows: 

1.  Compute  film  thickness  distribution  of  the  pad  which  is 
strictly  a  function  of  the  pivot  clearance  and  pitch  angle. 

2.  Integrate  pressure  to  obtain  pitch  torque  on  the  pad  about 
the  pivot  position. 

3.  Integrate  the  pad  dynamics  equations  through  one  time 
interval  to  obtain  the  pitch  of  the  pad. 

4.  Solve  Reynolds  equation  for  new  pressure  distribution. 

5.  Repeat  steps  1  through  4  until  the  pad  has  steadied  and 
pitch  angles  and  pressures  do  not  vary  through  success¬ 
ive  time  intervals. 
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6.  Print  load  capacity,  pivot  film  thickness,  minimum  film 
thickness,  pitch  angle,  friction  coefficient. 

By  selecting  proper  values  of  pad  inertia  this  process  can  converge  very 
rapidly.  Since  the  pivot  position  is  known  beforehand,  this  method 
obviates  the  necessity  of  producing  field  maps  and  obtains  steady-state 
data  very  expeditiously.  Using  this  technique,  single  pad  performance 
information  can  be  determined  as  a  function  of  pivot  clearance  for  a 
particular  set  of  conditions,  (A,  R/L,  <j>,  a,$). 

From  the  single  pad  information  full  bearing  performance  (three  pads) 
can  be  readily  computed.  Figure  33  shows  the  bearing  geometry.  For  the 
load  direction  between  pads , 

h  =  h  „ 
pi  P2 


h  =  G' 

Pi 


e  cosg 


=  =  r  [!_ 

n  n  *-  -L- 


e  cost 


Hp3  ^  [l  +  e']  for  no  spring 


H  _  =  £'  for  a  soft  spring  (stiffness  of 
^  C  spring  << 


stiffness  of  [2b] 

spring  << 

stiffness  of  bearing  film 


Using  the  single  pad  load  coefficient  vs  pivot  clearance  curve  and 
eq.  [1]  and  [2]  the  bearing  load  coefficient  can  be  determined. 

CL1  “  CL2 


CLT  "  2  CL1  “S3  "  CL3 


Note  for  the  spring  load  pad  is  also  obtained  from  the  single  pad 
information  but  it  must  also  equal  the  spring  force.  Therefore,  the 
spring  is  designed  using  as  its  load. 
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I 


JOURNAL  BEARING 


Figure  33.  Schematic  Representation  of  Pivoted-Pad  Journal  Bearing 


66 


For  the  load  direction  into  pad  number  3 

V  ’  V 

h  ,  -  C'  +  &'  cos8 

pl 

H  .  =  h  , 

Pi  _£i  =  —  [1  +  e"  cos(5] 

C  C 

Hp3  =  —  for  spring 

Hp3  =  -  [1  -  t']  for  no  spring 

c 

Therefore 

CLT  "  °U  '  2  CL1  C08B 
Total  friction  coefficient  (Fj)  is 


FT  =  F1  +  F2  +  F3 


The  minimum  clearance  for  the  bearing  is 


^  =  ^1  “  ^2 
^  =  ^3 


for  load  between  pads 
for  load  direction  at  pad  3 


[43 

[5a] 

[5b] 

[6] 

[7] 


Therefore  bearing  performance  can 
formation. 


be  obtained  from  the  single  pad  in- 


3.  SINGLE  PAD  PERFORMANCE  vs.  PIVOT  CLEARANCE 


68 


“  1.6 

o 


2  1.4 

Ui 


u3  1.2 

8 

<  1.0 

o 

_j 


A  *3.5 
R/L  -  0.25 


A  =  1.5 
R/L  =  1.0 


A=  3.5 
R/L  =  0.5 


A  =  1.5 
R/L  =  0.25 


S  0.8 


A=  1.5 
R/L  =  0.5 


A  =  3.5 
R/L  =  1.0 


PIVOT  FILM  THICKNESS,  Hp  -  hp/C 


Figure  34.  Single  Pad  Load  Coefficient  as  a  Function  of  Pivot  Clear- 
.  ance,  A  =  1.5,  3.5 


2.7 


Figure  36.  Single  Pad  Load  Coefficient  as  a  Function  of  Pivot  Clear 
ance,  A  =  10 
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PAD  LOAD  COEFFICIENT,  CT  -  W/p  RL 


-A  «20„  R,L*  .25 
■  A  =  20.,  R,  L  =  .5 
A  =  20.,  RyL  *  1.0 


PIVOT  FILM  THICKNESS,  H  »  h.  /C 
P.  P 


Figure  37.  Single  Pad  Load  Coefficient  as  a  Function  of  Pivot  Clear 
ance,  A  =  20 


FRICTION  COEFFICIENT 


Single  Pad  Minimum  Clearance  as  a  Function  of  Pivot 
Clearance,  A  =  1 .5,  3.5 


.  i  iiiimiiijj.,  ij, ,  jpp.«5* iLPt, .  111  iij  1 1  ■ij.Hiiiimnp’ 


o.o  0.2  0.4  0.6  0.8  1.0  1.2  1.4 


PIVOT  FILM  THICKNESS,  H  =  h  /C 
P  P 


Figure  43.  Single  Pad  Minimum  Clearance  as  a  Function  of  Pivot.  Clear 
ance,  A  =  20 
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4.  BEARING  LOAD  COEFFICIENT  AND  PIVOT  FILM  THICKNESS 
vs.  ECCENTRICITY  RATIO 
(ALL  PADS  RIGID) 
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BEARING  LOAD  COEFFICIENT,  CLT  -  W/P0RL 


1.4 


1.4 


Figure  44.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  1.5,  R/L  =  0.25 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  C.  T  -  W/P.RL 


.1 


Figure  45.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  1.5,  R/L  =  0.5 


PIVOT  FILM  THICKNESS,  Hp  =  hp/C 


BEARING  LOAD  COEFFICIENT  ,CLT-w/paRL 


1.4 


0.8  . 


A  •  1.5 

a  « ioo* 

R/L  ■  1.0 
fi  ■  60* 
<f>/a  -  0.65 


BEARING  ECCENTRICITY  RATIO,  E'-e/c' 

Figure  46.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Riqidlv 
Supported,  A  =  i.5,  R/L  =  1.0 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT 


c 

^ 

°',rt=^“7' 


J<v 

*Q£X 

c*o> 

Af  .  *» 


AF  <&' 


'/*  K/\ 
&°  <&?■ 


•s/<# 


®  %• 
/  G' 


A -3.5 

/  a «  ioo* 

R/L  -  0.25 
/9  -60* 

|  |  <£/a «  0.65  j 

0.2  0.4  0.6  0.8 

BEARING  ECCENTRICITY  RATIO,  E'-e'/c' 


Figure  47.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  3.5,  R/L  =  0.25 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT 


Figure  48.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  3.5,  R/L  *  0.5 
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PIVOT  FILM  THICKNESS 


BEARING  ECCENTRICITY  RATIO,  €'«  e'/c' 

wit 


a* 


Figure  49.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  3.5,  R/L  =  1.0 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  C._  -  W/P  RL 


✓  ✓ 

✓  / 

’  /  / 
S.  /  ' 


>*  *  ^  / 


Vl«  '/ 


V  *$$ 


A  -  5.0 
a  -  100° 
R/L  -  .25 

g  -  60° 

<t>/a  »  0.65 


BEARING  ECCENTRICITY  RATIO,  e'  -  e'/c ' 


. . H  - - 

- — 

c'/c»0.3 


Figure  50.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  5.0,  R/L  =  0.25 
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PIVOT  FILM  THICKNESS,  H  -  h  /C 


COEFFICIENT,  Ctt  =  W/P_ 


PIVOT  FILM  THICKNESS 


COEFFICIENT 


Figure  52.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  5.0,  R/L  =1.0 


I 

£ 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,-  CTT  «=  W/P  RL 


PIVOT  FILM  THICKNESS 


BEARING  ECCENTRICITY  RATIO,  e'  « 


Figure  54 


.  Bearing  Load  Coefficient  and  Pivot  Cl 
Ratio  for  Load  Direction  Between  Pads 
Supported,  A  =  10.0,  R/L  =  0.5 


BEARING  LOAD  COEFFICIENT,  C  -  W/P  RL 


PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT 


BEARING  ECCENTRICITY  RATIO,  e'  -  e'/c' 


Figure  56.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  20.0,  R/L  =0.25 
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BEARING  LOAD  COEFFICIENT 


Figure  57.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  20.0,  R/L  =  0.5 


1 5 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT 


Figure  58.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Direction  Between  Pads  with  All  Pads  Rigidly 
Supported,  A  =  20.0,  R/L  =1.0 
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PIVOT  FILM  THICKNESS,  H  =  h  /C 


Figure  59.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Supported 
A  =  1.5,  R/L  =  0.25 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  C^T  —  W/PaRL 


Figure  60.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Supported, 
A  -  1.5,  R/L  =  0.5 
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PIVOT  FILM  THICKNESS,  Hp  =  hp/C 


BEARING  LOAD  COEFFICIENT,  CLT  »  W/PqRL 


Figure  61. 


Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Supported, 
A  =  1.5,  R/L  =  1.0 


PIVOT  FILM  THICKNESS, 


BEARING  LOAD  COEFFICIENT,  CLT  -  w/paRL 


PIVOT  FILM  THICKNESS 


tricity 


BEARING  LOAD  COEFFICIENT,  C7  „  -  W/P  RL 


PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  C, _  -  W/P  RL 


LOAD  DIRECTED  AT  SHOE 


Figure  66.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Supported, 
A  =  5.0,  R/L  =  0.5 
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PIVOT  FILM  THICKNESS,  H  «  h  /C 


LOAD  DIRECTION  AT  PIVOT 


li  =  60° 
4>/c  =  0.65 


n 


\tjgZi 


BEARING  ECCENTRICITY  RATIO,  e'  «  e'  / c' 

Figure  67.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Supported 
A  ~  5.0,  R/L  ~  1.0 


PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  C  »  W/P.RL 


LOAD  DIRECTED  AT  PIVOT 


0.0  0.2  0.4  0.6  0.8  1.0 


BEARING  ECCENTRICITY  RATIO,  e'  =  e'/c' 


Figure  68.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Supported, 
A  =  10.0,  R/L  =  0.25 


PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  C.  »  W/P  RL 


LOAD  DIRECTED  AT  PIVOT 


A 


BEARING  ECCENTRICITY  RAIIO,  e"  -  e*/ c' 


Figure  70.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Supported, 
A  =  10.0,  R/L  =  1.0 
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PIVOT  FILM  THICKNESS,  H  -  h  /C 


BEARING  ECCENTRICITY  RATIO,’ z'  =  e'/c' 


WARING  LOAt)  COm'ICIKNT,  c,  •  W/P 


BEARING  ECCENTRICITY  RATIO,  c'  =  e'/c' 


Figure  72.  Beaming  Load  Coefficient  and  Pivot  Clearance  vs.  Ecce«-t 
Rat to  for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Si 
A  =  20.0,  R/L  =0.5 


108 


BEARING  ECCENTRICITY  RATIO,  e'  =  e'/c' 


Figure  73.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio^for  Load  Directed  at  Pad  with  All  Pads  Rigidly  Supported 
A  =  20.0,  R/L  =  1.0 


BEARING  LOAD  COEFFICIENT  AND  PIVOT  FILM  THICKNESS 
vs.  ECCENTRICITY  RATIO 
(ONE  SPRING  LOADED  PAD) 


BEARING  LOAD  COEFFICIENT 


Figure  74.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  1.5,  R/L  =  0.25 


PIVOT  FILM  THICKNESS 


Figure  75.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs*  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  1.5,  R/L  =  0.5 


BEARING  LOAD  COEFFICIENT 


0  0.2  0.4  0.6  0.8  1.0  1.2 

BEARING  ECCENTRICITY  RATIO,  €'  *  e'/c' 

Figure  76.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  «  1.5,  R/L  *  1.0 


U3 


THICKNESS 


BEARING  LOAD  COEFFICIENT,  C,T  *  W/p0RL 


/3  =  60* 
4>/a  *  0.65 
Ho-  *  c'/c 


C?/<o/  /  / 

'/  /  o /tofu  . 

8/  7  77°/ 


0.7 


0.6  * 


0.5  § 
o 


4e»  O' 


— |0.4  S 


0.2  0.4  0.6  0.8  1.0  1.2 

BEARING  ECCENTRICITY  RATIO,  €'  *  e'/c' 


Figure  77.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  3.5,  R/L  =  0.25 
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BEARING  LOAD  COEFFICIENT,  CLT  =  W/p0RL 


0  0.2  0.4  0.6  0.8  1.0  1.2 

BEARING  ECCENTRICITY  RATIO,  €'  =  e'/c' 

Figure  78.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Loadedf°A  !:°3d5DlJ®^ted0B®tween  Pads  with  Top  Pad  Spring 
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PIVOT  FILM  THICKNESS,  HD  =  hn/c 


a  *  100 


r/l  *  ;.o 


/3  -  60 


<p/a  =  0.65 


V  <0 
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UJ 
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BEARING  ECCENTRICITY  RATIO,  €'  *  e'/c' 


Figure  79.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  3.5,  R/L  =  1.0 


PIVOT  FILM  THICKNESS,  Hp  =  hp/c 


BEARING  LOAD  COEFFICIENT 


0  0.2  0.4  0.6  0.8  1.0  1.2 


BEARING  ECCENTRICITY  RATIO,  €'  =  e'/c' 


Figure  80.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  5.0,  R/L  =  0.25 
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PIVOT  FILM  THICKNESS,  Hp  =  hp/< 


BEARING  LOAD  COEFFICIENT,  CLT  =  W/p0RL 


Figure  81.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  5.0,  R/L  =  0.5 
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PIVOT  FILM  THICKNESS,  Hp  «  hp/c 


BEARING  LOAD  COEFFICIENT,  CLT  *  W/paRL 


Figure  82.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Spring 
Loaded,  A  =  5.0,  R/L  =  1.0 


119 


PIVOT  FILM  THICKNESS,  Hp  =  hp/< 


Figure  83.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  10.0,  R/L  =  0.25 
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BEARING  LOAD  COEFFICIENT,  CLT  =  W/pQRL 


il  Figure  84.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
*,*  Loaded,  A  =  10.0,  R/L  =  0.5 
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BEARING  LOAD  COEFFICIENT,  C 


0  0.2  0.4  0.6  0.8  1.0  1.2 

BEARING  ECCENTRICITY  RATIO,  €'  =  e'/c' 


Figure  85.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  10,  R/L  =  1.0 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  C,T  *  W/p0RL 


Figure  86.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for’  Load  Directed  Between  Pads  with  Top  Pad  Spring 
Loaded,  A  =  20,0,  R/L  =  0.25 
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PIVOT  FILM  THICKNESS,  Hp  =  hp/< 


BEARING  LOAD  COEFFICIENT,  C 


PIVOT  FILM  THICKNESS 


U  Figure  88.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  Between  Pads  with  Top  Pad  Spring 
; r  Loaded,  A  =  20,  R/L  =  1.0 
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PIVOT  FILM  THICKNESS, 


Figure  89.  Bearing  Load  Coeffii 
Ratio  for  Load  Dire' 


BEARING  LOAD  COEFFICIENT,  CLT  *  W/paRL 


Figure  90,  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  1.5,  R/L  =  0.5 
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FILM  THICK 


BEARING  LOAD  COEFFICIENT,  C,T  *  W/paRL 


8  5 

I 


Figure  91.  Bearing  Load  Coefficient  and  Pivot  Clearance  v$.  Eccentricity 

Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  1,5,  R/L  *  1.0 
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PIVOT  FILM  THICKNESS,  Hp  =  hp/c 


Figure  92.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  a  =  3.5,  R/L  =  0.25 


129 


PIVOT  FILM  THICK 


BEARING  LOAD  COEFFICIENT,  ClT  *  W/pcRL 


Figure  93.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  3.5,  R/L  =  0.5 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  CLT  =  W/pQRL 


Figure  94.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  3.5,  R/L  =  1.0 
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PIVOT  FILM  THICKNESS,  Hp  =  hp/c 


BEARING  LOAD  COEFFICIENT 


Figure  95.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  5.0,  R/L  =  0.25 
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PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT 


A  =  5.0 
a  =  100* 
R/L  =  0.5 

P  =  60° 
4>/a  =  0.65 
Hp  *  c'/c 
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BEARING  ECCENTRICITY  RATIO,  €'  =  e'/c' 


Figure  96.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  5.0,  R/L  =  0.5 


PIVOT  FILM  THICKNESS 


BEARING  LOAD  COEFFICIENT,  C,T  *  W/p0RL 


Figure  97.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 

Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  5.0,  R/L  =  1.0 
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PIVOT  FILM  THICKNESS,  Hp  =  hp/c 


PIVOT  FILM  THICKNESS, 


BEARING  LOAD  COEFFICIENT 


Figure  99.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  10,  R/L  *  0.5 
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PIVOT  FILM  THICKNESS,  Hp  «  hp/c 


BEARING  ECCENTRICITY  RATIO,  €'  *  e'/c' 

Figure  100.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  a  =  10,  R/L  =  0 


BEARING  LOAD  COEFFICIENT,  C 


j3  =  60* 
<£/a  »  0.65 


Hd_  *  c'/c 
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BEARING  ECCENTRICITY  RATIO,  €'  *  e'/c' 


Figure  101.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  20,  R/L  =  0.25 
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PIVOT  FILM  THICKNESS, 


BEARING  LOAD  COEFFICIENT,  CLT  *  W/paRL 


Figure  102.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Rat10  for  Load  Directed  at  Spring  Loaded  Pad,  A  =  20,  R/L  =  0.5 

■m* 
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PIVOT  FILM  THICKNESS,  Hp  =  hp/< 


BEARING  LOAD  COEFFICIENT,  CLT  *  W/pQRL 


Figure  103.  Bearing  Load  Coefficient  and  Pivot  Clearance  vs.  Eccentricity 
Ratio  for  Load  Directed  at  Spring  Loaded  Pad,  A  »  20,  R/L  =  0.5 
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APPENDIX  II 

VARIABLE  GRID,  DU  FORT-FRANKEL  ANALYSIS  OF 
GAS-LUBRICATED  BEARINGS 

1.  INTRODUCTION 

Solution  of  high  A  gas-bearing  problems  have  posed  numerical  problems 
because  of  the  very  high  pressure  gradients  at  the  regions  of  minimum 
film  thickness.  Use  of  variable  grid  methods  enable  a  dense  grid  at 
regions  of  high  gradients,  and  a  coarse  grid  otherwise  and  circumvents 
the  problems  of  numerical  instabilities  imposed  by  uniform  grids  of 
spacing  that  would  not  absorb  excessive  computer  storage  or  computational 
time. 

Another  problem  inherent  in  time  transient  lubrication  problems  is 
repetitive  solution  of  the  Reynolds'  equation  throughout  the  grid  at  each 
time  increment.  Implicit  formulations  are  encumbered  by  large  numbeis  of 
matrix  inversions  while  the  usual  explicit  formulation  requires  very 
small  time  increments  for  numerical  stability.  The  DuFort-Frankel  Method 
combines  the  advantages  of  explicit  solution  with  reasonable  time  incre¬ 
ments.  Thus,  the  discussion  following  incorporates  two  important 
features  with  regard  to  gas  bearing  analysis,  namely 

a)  Variable  grid 

b)  Stable  explicit  solution  for  pressure  with  practical  time 
increments . 


2.  ANALYSIS 

Making  all  the  standard  assumptions  for  compressible  lubrication 
theory  (contact  viscosity,  negligible  inertia  effects,  perfect  gas,  etc.) 
the  continuity  eq  can  be  written 

4  /  pdv  +  /  p  V  •  n  ds  =  0  [1] 

3t  v  K  s  y  1  J 

where 

V  "  V0  \  +  Vz  nz 
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_  1  3p  h  BO 

R  30  12y  2 


v  .-iEliL- 

VZ  3Z  12  y 

dv  =  hRd0dz 


By  non-dimensionalizing  [13  with 


HC  -  h 
ZL  -  z 


and  letting  $  »  constant 


,  1  3p_  P3p 

and  vr  t~~  ~ 

2  3n  3n 


equation  [1]  becomes 


3  2 

3  IS  HPd0dZ  -  I  [(?-  |£~  "  A  pH)dz 
3T  2  39 

2  3  2 

,  HJ  3T 

+  (L>  T  3Z~  d0] 


Integrations  are  accomplished  by  accumulation  of  the  net  flow 

emanating  from  a  closed  region  surrounding  each  grid  point.  The  equation 

for  each  grid  point  is  made  algebraic  by  using  finite  increments  for 

distances  and  two  point  finite  differences  for  derivatives.  The  region 

for  variable  grid,  shown  dotted  on  Figure  104,  encompasses  half  the  dis- 

th 

tance  between  the  i,  j  point  and  its  neighbors. 


The  Du-Fort-Frankel  technique  is  introduced  by  substituting  for 

t*V» 

the  pressure  at  the  i,  j  point  at  time  t,  the  average  of  future  and 
past  pressures,  i.e. : 


The  formulation  of  the  finite-difference  continuity  equation  is  a 
quadratic  of  the  form: 


+  *  t  +  1 

APij  +BPij  +c  =  ° 


where  the  coefficients  are  functions  of  the  local  film  thickness  and 
grid  spacing,  compressibility  parameter  and  time  increment.  It  can  be 
shown  that  as  At  •+  0  the  pressure  solution  to  equation  [4]  converges  to 
the  correct  solution.  Solution  of  the  quadratic  equation  [4]  at  each 
grid  point  is  numerically  stable  for  a  reasonably  large  time  increment 
and  computational  time  is  much  less  than  either  the  usual  implicit  or 
explicit  schemes  described  by  References  4  and  5. 
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APPENDIX  III 

DYNAMIC  ANALYSIS  AND  DESIGN  CURVES 
FOR  TILTING  PAD  JOURNAL  BEARING 


1.  STABILITY  MAPS 


SHOE  PITCH  INERTIA  PARAMETER  I 


*-*jv  .'~*v  «<*.-  >~<  ‘“Kw^i'^a-c^®J£^Lf 


Figure  105.  Stability  Map,  R/L  =  0.5,  A  =  1.5 
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SHOE  PITCH  INERTIA  PARAMETER  I 


Figure  106.  Stability  Map,  R/L  =  1.0,  A  =  1.5 


a  =  100° 
R/L  =  0.5 

<f>/a  -  0.65 
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Figure  107.  Stability  Map,  R/L  =  0.5,  A  =  3.5 
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STABLE 


A  «  3.5 
a  *  100° 
R/L  *  !.0 
<f>/a  »  0.65 


J L  1  -1-1  llll _ I LJJJLLiil _ L-J_LJ 

10  too 

SHAFT  MASS  PARAMETER  Ms  -  RL 


Figure  108.  Stability  Map,  R/L  =  1.0,  A  =  3.5 
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Figure  109.  Stability  Map,  R/L  =  0.5,  A  =  5.0 
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SHOE  PITCH  INERTIA  PARAMETER 


Figure  110.  Stability  Map,  R/L  =  1.0,  A  =  5.0 
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SHOE  PITCH  INERTIA  PARAMETER  I 


Figure  111.  Stability  Map,  R/L  =  0.5,  A  =  10 
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SHOE  PITCH  INERTIA  PARAMETER  I 


Figure  112.  Stability  Map,  R/L  =  1.0,  A  =  10 
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SHOE  PITCH  INERTIA  PARAMETER  I 


Figure  113.  Stability  Map,  R/L  =  0.5,  A  -  20 
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SHOE  PITCH  INERTIA  PARAMETER  I 


Figure  114.  Stability  Map,  R/L  =  1.0,  A  =  20 
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